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ABSTRACT 



We use the algorithm of ICole et alJ 1)20001) to generate merger trees for the first 
star clusters in a ACDM cosmology under an isotropic UV background radiation field, 
parametrized by J21. We have investigated the problem in two ways: a global radiation 
background and local radiative feedback surrounding the first star clusters. 

Cooling in the first halos at high redshift is dominated by molecular hydrogen, 
H2 — we call these Generation 1 objects. At lower redshift and higher virial tempera- 
ture, T vir > 10 4 K, electron cooling dominates — we call these Generation 2. 

Radiation fields act to photo-dissociate H2, but also generate free electrons that 
can help to catalyse its production. At modest radiation levels, J21 / (l + z) 3 ~ 10~ 12 — 
10~ 7 , the nett effect is to enhance the formation of Generation 1 star-clusters. At 
higher fluxes the heating from photo-ionisation dominates and halts their production. 
With a realistic build-up of flux over time, the period of enhanced H2 cooling is so 
fleeting as to be barely discernable and the nett effect is to move primordial star cluster 
formation towards Generation 2 objects at lower redshift. 

A similar effect is seen with local feedback. Provided that enough photons are 
produced to maintain ionization of their host halo, they will suppress the cooling in 
Generation 1 halos and boost the numbers of primordial star clusters in Generation 2 
halos. Significant suppression of Generation 1 halos occurs for specific photon fluxes 
m excess of about lO^phs^Mo -1 . 
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1 INTRODUCTION 

Primordial star clusters contain the first stars to form 
in t he Universe, from zero-metallicity gas. Previous work 
( e.g.lTegmark et alJll997t IXb el et a flil998tlHutchings et alJ 



well have a negative feedback effect on subsequent popula- 
tion III star formation by supressing cooling via this mecha- 



12002 - iBromrn et aljl999l:lsantoro fc Thomas 2003) has con- 
centrated on the very first objects for which there is no sig- 
nificant external radiation field. However, these first clusters 
are expected to produce massive stars which will irradiate 
the surrounding Universe and may well be responsible for 
partial re-ionisation of the intergalactic medium. 

The only species produced in sufficient abundance to 
affect the cooling at early times is molecular hydrogen. Its 
presence allows the first objects to cool and form in low tern- 



nism. This problem is not a trivial one and has been the sub- 
ject of much interest in recent years (e.g. Haiman ct al. 1996; 


Haiman et al J 12000: 


Kitavama et al. 2000; Glover & Brand 


2001^ 


lOmukal 12001 


iMachacek et al. 


120011: iRicotti et alJ 


2002 


Oh & Haiman 


2002: ICenl 12003 


: ICiardi et alJ 12003: 


Yoshida et all 12003. 


Tumlinson et alJl2004b The complex- 



perature halos (T < 10 K) at high redshift (2 



20-30). 



However, molecular hydrogen is very fragile and can easily 
be dissociated by UV photons in the Lyman- Werner bands, 
(11.2 — 13.6eV). Thus, the formation of the first stars may 
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ity of the feedback and the large number of unknowns (e.g. 
population III IMF, total ionising photon production, etc.) 
make this problem very challenging. 

In an attempt to understand this era of primordial 
star cluster formation, we in vestigated in a previous pa- 
per iSantoro fc Thoma d 120031 . hereafter ST03) the merger 
history of primordial haloes in the ACDM cosmology. There 
we assumed no external radiation field (other than that pro- 
vided by CMB photons). The Block Model of lCole fc Kaiserl 
( 1988) was used to generate the merger history of star clus- 
ters using a simple model for the collapse and cooling crite- 
rion, hence identifying those haloes that were able to form 
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stars before being disrupted by mergers. We then contrasted 
the mass functions of all the resulting star clusters and those 
of primordial composition, i.e. star clusters that have not 
been contaminated by sub-clusters inside them. We found 
two generations of primordial haloes: low-temperature clus- 
ters that cool via H2, and high-temperature clusters that 
cool via electronic transitions. 

We investigated two regions of space each enclosing a 
mass of 10 11 /i _1 Mq: a high-density region corresponding to 
a 3<r fluctuation (So = 10.98), and a mean-density region 
(So = 0), where So is the initial overdensity of the root block. 
In the high-density region we found that approximately half 
of the star clusters are primordial. The fractional mass con- 
tained in the two generations was 0.109 in low temperature 
clusters and 0.049 in high temperature clusters. About 16 
per cent of all baryons in this region of space were once part 
of a primordial star cluster. In the low-density case the frac- 
tional mass in the two generations was almost unchanged, 
but the haloes collapsed at much lower redshifts and the 
mass function was shifted toward higher masses. 

In the present paper, as a continuation of the previous 
work, we include the effect of ionising radiation in two dif- 
ferent ways: firstly we add a homogeneous background radi- 
ation field; secondly, we consider feedback from the first star 
clusters formed in the merger tree — these will form an ionis- 
ing (and photo-dissociating) sphere around them, changing 
the cooling properties of neighbouring star clusters. A fur- 
ther improvement upon our previous work includes the use 
of a more realistic merger tree. 

We describe our chemical network including radiative 
processes in Section [5] and the new merger tree method 
in Section |3] The effect of a global ionisation field on the 
formation of stars in primordial star-clusters is considered 
in Section 2] and that of local feedback in Section [S] Finally, 
we summarise our conclusions in Section |H] 



2 PRIMORDIAL CHEMISTRY AND GAS 
COOLING IN THE PRESENCE OF 
RADIATION 

2.1 Chemical model 

In this Section we introduce the chemical network needed to 
follow the coupled chemical and thermal evolution under a 
homogeneous UV background radiation field. 

The non-equilibrium che mistry code is ba s ed up on the 
minimal model presented in iHutchines et all £2002, here- 
after HSTC02). It calculates the evolution of the follow- 
ing 9 species: H2, H, H + , Hj, H~, He, He + , He ++ and e~. 
The important cooling processes are: molecular hydrogen 
cooling, collisional excitation and ionisation of atomic Hy- 
drogen, collisional excitation of He + , and inverse Compton 
cooling from cosmic microwave background photons. In this 
paper we only consider the low density-high temperature 
(T > 300K) limit. Thus, we have ignored the effects of HD 
cooling, which is only important in the high density-low tem- 
perature regime. 

To this chemical model we have added the photo- 
ionisation and p hoto-dissociation reactions compile by 
lAbel et alJ jl997f) and listed in Table □ This consists of 9 
reactions involving the interaction of each species with the 



background photons. These are: photo-ionisation of H, He, 
He + and H2, with threshold energies of hv =13.6, 24.6, 54.4 
and 15.42 eV, respectively; photo-detachment of H~with a 
threshold energy of 0.755 eV, potentially an important pro- 
cess since H~catalyses the formation of H2; and photo- 
dissociation of H^" and H2 (by the Solomon process and 
by direct photo-dissociation). In the case of the Solomon 
process dissociation happens in a very narrow energy range 
12.24eV< hv < 13.51eV. 

The energy equation takes the following form: 



d(n t T) 



(A he 



Ac 



(1) 



dt 3k 

where nt is the total number density of all species, T is the 
temperature, A coo i and Aheat are the cooling and heating 
terms, respectively, and we have assumed a monatomic en- 
ergy budget of |fcT per particle (the energy associated with 
rotational and vibrational states of H2 is negligible). 



A coo l = Ah.cc + AH.ci + A H<J + ce + Ah 2 ,cc + Aco 



(2) 



where the suffixes ce and ci mean collisional excitation and 
ionisation respectively, and expressions for each of these 
terms are given in HSTC02. 



Aheat — A 



H,pi + Anc.pi + A Hc + pi + AH 2 ,pi + AH 2 ,pd 



(3) 



where the suffixes pi and pd mean heating from photo- 
ionisation and photo-dissociation respectively. The expres- 
sion for each term was calculated using Equation IB1I from 
the Appendix, using the cross-sections listed in Tabled 



2.2 UVB spectrum 

The non-equilibrium chemistry and the thermal evolution 
of the clouds are calculated assuming the presence of an 
ultraviolet radiation field of the power-law form 



1r ,-21 -1 -2tt -1 -1 

10 erg s cm Hz sr , 



(4) 



where hvn = 13.6 eV is the Lyman limit of H. Here the 
direction of the radiation field is not important and the nor- 
malisation is given in terms of the equivalent isotropic field. 

In the present-day Universe, QSOs hav e much steeper 
spectra (a « 1.8, e.g. |Zh^ri£_et_alJ 19971) than do stars 
(a w 5, e.g. iBarkana^^Loe bl 119991) . However, Popula- 
tion ill stars are likely to be biased to much higher spec- 
tral energies and their spectra may resemble those of 
QSOs above the Lyman limit jTumlinson fc Shu J l2000l : 
iTumlinson. Shull fc Venkatesanl l2003fi . In this paper, we 
take a = 2 which could equally well apply to either type 
of source. 

For the calculation of all the photo-ionisation and 
photo-dissociation rates and as well as the heating terms, 
we assume an optical thin medium and no self-shielding. 



2.3 Cooling of isolated clouds 

Apart from a scaling factor, the cooling time, t C ool, of iso- 
lated halos subject to a uniform radiation field depends 
solely on the ratio of the number densities of baryons, rib, 
and photons, n 7 , i.e. t coo i =fn(T v i r ,nb/n~ f )/nb- Thus the 
effect of the global background radiation on single haloes 
can be presented in two ways: we can fix either the baryon 
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Table 1. This table summarises the important reactions that should be included in a chemical network if a uniform background radiation 
field is present. Compiled by Abel et al. (1997), except cross sections 25 and 27. The number index of each reaction corresponds to those 
in that paper. Reference: Osterbrock (1989, 089), de Jong (1972, DeJ72), O'Neil & Reinhardt (1978, OR78), Tegmark et al. (1997, 
TSR97), Haiman, Rees & Loeb (1996, HRL96), Abel et al. (1997, AAZN97). 



Reaction 



cross sections/cm 2 



Reference 



20 H + 7 i 



21 He- 



26 
27 
28 



7 i 



H+ + 2c~ 
-+ He+ + e 



22 He+ +71 — > He++ - 

23 H-+71 >H + c- 



24 H 2 



25 IE 



-71 



H 2 +7 ' 
H 2 +7 >- 
H 2 +7 h 



H, 



» H + H+ 

» 2H+ + e 
H 2 * 1— 
H + H 



f21 
o"22 
C23 



C"25 



7.42 x 10~ 18 



A = 6.30 x 10 _18 cm 2 



(l_e-=™/e) J I e =y / u /v th - 1, hi/ tA = 13.6eV 

/ / „ \ -2.05 / „ \ -3.05\ 

(A /^ 2 )(^ : ) _4 ( e( ^"w'r )' /,,/th = z2 x i3,6ey and z = z 

7.928 x 10 5 (i/ - v th )i (4;-) , hi/ > hv th = 0.755eV 





6.2 x 10- 18 hu - 9.40 x lO" 17 
1.4 x 10~ w hv - 1.48 x 10~ 17 
2.5 x 10- 14 (^)- 2 ' 71 

7.401 x io- ls lo( _:r2_ 0302a:3_ ' 0158:l:4 ) 



hv < 15.42eV 
15.42 <hu < 16.50eV 
16.5 <hv < 17.7eV 
hv > 17.7eV 



2.762 ln(Ai//11.05eV) 
hv > 2.65eV 



1Q-16. 926-4.528x10 A hi/+2.238 x 10 » (hi/) z +4.245x 10 30eV< /ll^ <90eV 



H + H cr 27 = 3.71 x 10- 



12.24eV< hf <13.51cV 



see Reference for the expression 
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Figure 1. The cooling time, i coo i in the J21 — T v i r plane. The 
colour bar is logl0(t coo i/yr). This plot was calculated at a fixed 
density corresponding to z = 20. 
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Figure 2. The molecular hydrogen fraction, /h 2 ; after one cool- 
ing time, in the J%\ — T v i r plane. The region corresponding to 
clouds that have not cooled by the time the expansion factor 
reaches 10 has been coloured white. 



density or the amplitude of the radiation field and vary the 
other. 

Figure shows the cooling time, t coo i, of haloes ex- 
posed to different levels of background radiation for a fixed 
density of nn ~ 0.31cm -3 , corresponding to the mean 
density within a collapsed halo at z = 20, as seen from 
the J21 — Tvir plane. t coo i is defined as the time that the 
halo takes to cool from virialisation to the moment when 
T = 0.75T v ir (or, if it doesn't cool, we stop the integration 
when the scale factor reaches 10) . Any halo that falls on the 
red region will not be able to cool in a Hubble time. On the 



other hand haloes on the dark blue or violet part of the plot 
will cool in much less than a dynamical time. 

In Figure |21 we show the fractional density of molecular 
hydrogen, /h 2 , after one cooling time in the J21 — T v i r plane. 
The region that has not cooled by the time the scale factor 
reaches 10 has been coloured white. For the zero-flux case, 
we see that the highest molecular fractions occur for a virial 
temperature of just over 10 4 K. Halos of this temperature are 
able to partially ionize hydrogen atoms and the free electrons 
then go on to catalyse production of H2. 

From these plots, we can see that the ionising radiation 
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has two main effects. Firstly, it provides a heat source that 
prevents cooling of halos. The minimum virial temperature 
of halos that can cool gradually increases along with the 
photon flux, whereas the heating has little effect on halos 
above this temperature. The effect is visible in Figure as 
a sharp transition between the green/blue and red regions 
that runs diagonally across the plot. 

Secondly, the ionizing flux can boost cooling of low- 
temperature halos (Tvi r < f 4 K) by creating free electrons 
that then catalyse the production of H2. For any given tem- 
perature, there is only a narrow range of photon fluxes for 
which this is important before photo-heating swamps the 
increased H2 cooling. 

Thus, as the background ionising radiation builds up in 
the Universe, we expect there at first to be a small boost 
in the formation of low-temperature (low-mass) star clus- 
ters and then a sharp decline. The precise values of J21 at 
which this occurs will depend upon the redshift of struc- 
ture formation (recall that the above plots are for a red- 
shift of 20 and the required values of J21 will scale as 
(f + z) 3 ). In Section 2] we consider both specific values of 
J21 = I0 -10 , 10 -5 , 10 _2 & 10 and also a time-evolving ioni- 
sation field. 



3 MERGER TREE 

In a previous paper JSantoro fc Thomasi20o1) we generated 
m erger histories o f dark matter halos using the block model 
of ICole fc Kaiser] <u988l) . In this model, a parent block of 
mass Mo and density fluctuation 80 is halved producing two 
daughter blocks of mass Mi — Mo/2. Extra power, drawn at 
random from a Gaussian distribution, is then added to one 
block and subtracted from the other in order to conserve the 
overall level of fluctuations in the root block. This process is 
then repeated until the desired resolution has been reached. 
The reader should note that the tree is produced by stepping 
back in time. Therefore, the formation of the daughter blocks 
occurs at a higher redshift than their parent block. A valid 
criticism of this model is that the mass distribution is not 
smooth, as each level changes in mass by a factor of two. 
This discretisation of mass is an undesirable constraint on 
our model and for this reason we have moved to a more 
realistic method of generating the merger tree. 

We now use a Monte-Carlo algorithm to gener- 
ate the merger tree. There are a number of codes 
which use this technique (e.g. iKauffmann fc White! Il993t 
ISomerville fc Kolattl I1999T) how e ver, w e have chosen the 
method proposed bv lCole et alJ (120001) . In short, the code 
uses the extended Press- Schecter formalism to generate the 
merger histories of dark matter halos, but we refer the reader 
to this paper for a complete discussion of the code. The main 
advantage of this method over the block model is that the 
discretisation of mass has been removed. 

Unlike most previous tree implementations, we do not 
place our merger tree onto a predefined grid of timesteps but 
utilise a very fine time resolution. This allows us to retain 
the block model's prescription of two progenitors per merger, 
consistent with our previous work. The algorithm allows us 
to account for the accretion of mass below the resolution 
limit, however we assume that this will not seriously affect 
the structure and cooling of the halo. 



An important aspect of the code is how we treat the 
mergers of halos. The outcome of a merger will depend upon 
the ratio of masses of the merging halos, 



where Mi < M2. Consequently, we introduce, </ m m, as the 
minimum mass ratio to affect the cooling of a halo. This 
leads to two cases: 

• If q > g m i n then the two halos merge and their cooling 
is completely disrupted. The gas is then shock heated to the 
virial temperature of the parent halo, erasing all previous 
cooling information, and the cooling starts afresh. 

• If q < Qmin then we assume that the smaller of the 
two daughters is disrupted. We then compare the times at 
which the larger daughter and the parent halo would cool. 
If the former occurs first then we postpone the merger and 
allow the cooling of the daughter to proceed; otherwise we 
continue as for q > qmin- 

In the future we would like to determine an appropri- 
ate value for q m in from hydr odynamical simulations. S imula- 
tions of galaxy mergers fe.g. lMihos fc Hernauistlf 996) show 
that when objects with mass ratios q > 0. collide the galactic 
structure of both objects is seriously disrupted. They clas- 
sify these as 'major' mergers. We suspect that smaller mass 
ratios would still sufficiently disrupt the cooling gas cloud. 
Consequently, for this work we set q m i n = 0.25. Although we 
do not show it here, our results remain largely unchanged in 
the range 0.2 < q m j n < 0.3. 

We treat the metal enrichment of halos in the same way 
as ST03 where it was assumed that, regardless of whether 
or not star clusters survive a merger, they instantly contam- 
inate their surroundings with metals and the enrichment is 
confined to the next level of the merger hierarchy (i.e. the 
metals do not propogate into halos on other branches) . This 
is the same for both the global (sectional and local models 
(section I^J. 

Once a halo has been contaminated it is no longer 
classed as primordial, irrespective of whether it can form 
more stars or not. However, no attempt is made to account 
for the transition from population III to population II star- 
formation. Consequently, contaminated halos are assumed 
to cool at the same rate as their primordial counterparts 
and in the case of our local model, produce the same ionis- 
ing flux. We intend to investigate the effect of this transition 
in future work by including cooling from metals. 

In this paper, in common with ST03, we use a root 
mass of IO n M for our tree, and a mass resolution of 
9.5 x 10 4 Mq. However, we use slightly differen t cosmological 
param eters as derived from the WMAP data llSpergel et alJ 
12003ft of n = 0.3, A = 0.7, Q b0 = 0.0457, h = 0.71, 
as = 0.9, and a power spectrum as calculated by CMBFAST. 

Figure |3] shows a comparison between the new merger 
tree (red line) and the older Block Model from ST03 (blue 
line) of the fractional mass per dex of primordial star clus- 
ters as a function of virial temperature, averaged over a large 
number of realisations and in the absence of ionising radia- 
tion. It is clear that the new tree has had a significant affect 
on the number of primordial objects that are formed. While 
the total number of objects that are able to cool remains 
roughly constant, our new model produces only a third of 
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Figure 3. The fractional mass per dex of primordial star clusters 
as a function of virial temperature: red line - new merger tree, 
blue line - Block Model. 



1 : 3 < z < 5 
(i±£)4 : < z <3. 

We adopt their model and fiducially fix the onset of the UV 
background at zuv = 50 (at which time the normalisation is 
negligible) . We have added a factor of /3 into this expression 
so as to control the rate at which the field builds up and 
present results for 3 cases: (3 = 0.8 (rapid), 1 (standard), 
and 2 (slow). 

For the latter cases especially, the global ionising flux 
can be thought of as coming from pre-existing star clusters 
(or quasars) that form in high-density regions of space and 
that are gradually ionising the Universe around them. For 
this reason, we take the mean-density of the tree to be equal 
to that of the background Universe. In Section we will 
consider a high-density region for which the ionisation field 
is generated internally from the star clusters that form in 
the tree. 



the primordial objects compared with the original, although 
the mass fraction has only been reduced by half. Qualita- 
tively the results remain unchanged: we still observe two 
generations of halos, distinguished by their primary cooling 
mechanism, as discussed in ST03. In addition, we have re- 
moved all features associated with the discrete mass steps 
(e.g. the feature at ~5000K in the original model). 

The smoothed mass distribution has lead to many un- 
equal mass mergers which were not present in the previous 
model thus increasing the likelihood of contaminating large 
haloes with much smaller ones which happen to cool first. 
Equally, the chance that haloes are involved in mergers that 
disrupt their cooling is increased. These effects conspire to 
reduce the overall number of primordial objects. 



4 GLOBAL IONISATION FIELD 
4.1 Model 

In this section, we consider the effect of a global ionisation 
field that affects all halos equally. As previously mentioned, 
we will restrict ourselves to a power law ionising flux with 
index a — 2, corresponding either to a quasar spectrum or 
that of stars of primordial composition. We present results 
for four cases of constant normalisation: J21 = 10~ 10 , 10~ 5 , 
10" 2 & 10. These are chosen to be representative of a very 
low flux where the effect on each halo is minimal; a flux 
which has a positive effect on the capacity of the gas to 
form H2; and two examples of higher amplitude fluxes that 
destroy H2. 

In addition to the cases outlined above, we consider a 
time- dependent build up of the background flux. Using ra- 
diative hydrodynamical calculations to examine the effect of 
the UV background o n the collapse of pre-galactic clouds, 
iKitavama et al.l <|2000T) modelled the evolution of a UV back- 
ground as: 



4.2 Results 

In Figure |U we plot the fractional mass of primordial star 
clusters as a function of (a) virial temperature and (b) halo 
mass. The colours red, blue, cyan & green correspond to 
J21 = 10" 10 , 10" 5 , 10~ 2 & 10, respectively— note that the 
peak of the distributions do not move steadily from left to 
right as J21 is increased. Figure shows histograms of the 
star-formation redshifts of the primordial clusters. 

The lowest amplitude case is almost indistinguishable 
from that of zero flux. For this reason we have not plot- 
ted the latter. There are two bumps in the virial tempera- 
ture histogram corresponding to two distinct cooling mech- 
anisms. In ST03, these were christened Generation 1 (low- 
virial temperature, T < 8 600 K, low-mass, high collapse 
redshift, dominated by H2 cooling) and Generation 2 (high- 
virial temperature, T > 8 600 K, high-mass, low collapse red- 
shift, dominated by electronic cooling). 

As the flux is increased, the effect of the radiation field 
is to promote the cooling of Generation 1 halos. The typ- 
ical virial temperature and mass of such halos decreases, 
and the number of Generation 2 star clusters is reduced as 
collapsing halos are more likely to have been polluted by 
metals from smaller objects within them. For J21 = 10 -5 
(blue curve), the effect is so pronounced that it completely 
eliminates Generation 2 objects. However this is an extreme 
case, because, as FiguresQ&I^Jshow, this flux has been cho- 
sen to produce close to the maximum possible enhancement 
in the H2 fraction and a corresponding reduction in cooling 
time throughout the redshifts at which these halos form. 

If the background flux is increased further, then the 
enhancement in Generation 1 star clusters is reversed. For 
J21 = 10 -2 (cyan curve) the balance has shifted almost en- 
tirely in the favour of Generation 2 clusters and by J21 = 10 
all Generation 1 clusters have been eradicated. 

It is interesting to note that the mass fraction of stars 
contained in primordial star clusters is not greatly affected 
by the normalisation of the ionising radiation, varying from 
0.05 to 0.1. However, the mass (and hence number) of the 
star clusters varies substantially. A modest flux will increase 
the number of small clusters, moving the mass function to 
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(c) 




Moss (M e ) 



Figure 4. Fractional mass per dex of primordial objects as a 
function of (a) virial temperature and (b) mass, for 4 different 
cases: red corresponds to J21 = 10 — 10 ; blue to J21 = 10 — 5 ; cyan 
to J21 = 10~ 2 ; and green to J21 = 10. 




5 10 15 



Figure 5. Histograms of star-formation redshifts for primordial 
halos. The colour coding is the same as in Figure HI 

virial-temperature star clusters, is more complicated. In this 
paper we are concerned with primordial objects, by which 
we mean those with zero metallicity. The key question, then, 
is whether large halos are contaminated by metals from sub- 
clusters within them. With a slow build-up of flux cooling in 
these sub-clusters is enhanced resulting in increased contam- 
ination and a reduction in the number density of primoridal 
Generation 2 halos. However, a rapid build-up of flux cuts 
off production of small halos dramatically and the number 
of primordial Generation 2 halos is increased. 



5 LOCAL FEEDBACK 

In this section, we consider local feedback, i.e. that produced 
by star clusters internal to the tree. 




5.1 Model 



lower masses while a greater flux produces the opposite ef- 
fect. 

The number of primordial star clusters as a function 
of star formation redshift is shown in figure The redshift 
distribution is similar in all cases, which differs significantly 
from model with time- varying flux, discussed next. 

Figure |S| plots the fractional mass of primordial star 
clusters, as a function of (a) virial temperature and (b) 
mass. Unlike our previous results, the introduction of a time- 
evolving field has had a devastating effect on the mass frac- 
tion of primordial objects. As the rate at which the flux 
builds-up increases, the peak of the Generation 1 mass func- 
tion moves towards lower masses (and virial temperatures) 
and its normalisation decreases significantly. At the same 
time, as shown in figure Q the peak in the production rate 
moves to higher redhsifts. In each case, it corresponds to a 
flux of J21 ~5x 10~ 5 for which, from the previous results, 
an enhancement in the H2 fraction and hence a dcrease in 
the cooling time is expected. 

The effect on the production of higher-mass, higher- 



5.1.1 Physical picture and assumptions 

The large value of the Thomson electron-scat tering optical 
depth in the WMAP data of r = 0.17 ± 0.04 jSpergel et alJ 
2003) suggests an early reionisation era at 11 < 2 < 
30 (180±go°Myr after the Big Bang). This requires a high 
efficiency of ionising photon production in the first stars, cor- 
responding to a deficit of low-mass stars. This i s backed up 
by theoretical arguments and simulations fe.g. lAb el et al* 



20021: iBromm et al]|2002l ISchaered 120021 : iTumlinson et al 
200i ISantoro fc Shullll2005l7 of star-formation in a primor- 



dial gas, for which fragmentation was found to be strongly 
inhibited by inefficient cooling at metallicities below about 
10- 3 ' 5 Z Q . 

There is still some debate as to whether the first 
star f orming halos wi ll produce a single massive star 
(e.g. lAbel et al] 120021) or fragment further to form the 
first star clusters (e.g. IBromm et al] Il999l) . Whichever 
of these is correct ma kes little difference to our results. 
ITumlinson et all i2004f) considered a number of IMFs that 
may lead to the required early reionisation. These have 
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Figure 6. Fractional mass per dex of primordial objects as a 
function of (a) virial temperature and (b) mass, for 4 different 
cases: orange corresponds to the no flux case; blue to fi = 2; red 
to /3 = 1; and green to /3 = 0.8. 



ionising fluxes per unit mass in the range Q « 10 - 
10 48 phs -1 Mq -1 , but they all have similar ionising efficien- 
cies of about 80 000 photons/baryon when integrated over 
the whole life of the stars. A value which is similar to that 
obtained for single massive stars [M > 2OM0 ) . If we assume 
a uniform flux over time, then this corresponds to mean life- 
times of 3.0 x 10 7 -3.0 x 10 6 yr, respectively. Note that these 
values of Q are much greater than t he average for the Milk y 
Way, Q » 8.75 x 10 43 phs" 1 Mq" 1 jRicotti fc Shuilll200CD . 

When halos in the merger tree are able to cool, we as- 
sume that they will convert part of their baryonic compo- 
nent into a "star cluster" (primordial or otherwise). These 
objects will exert radiative feedback onto the next genera- 
tion of haloes that form inside the same tree. The photon 
flux will also depend upon the star-formation efficiency and 
the escape fraction from the star-forming region in the cen- 
tre of the halo. In this paper, we are not concerned with the 




Figure 7. Histograms of star-formation redshifts for primordial 
halos. The colour coding is the same as in Figure IdI 



Table 2. Parameters of the ionisation models that we consider: 
model number; fraction of mass in stars, /„; specific ionising 
fluxes, Qoi i n units of phs -1 M0 -1 ; total number of ionising 
photons, Wo, and ionising flux, So, in units of phs -1 per solar 
mass of matter (baryonic plus dark); and the time for which the 
ionising flux acts, t*, in units of years. 



id 


/. 


Qo 


N 


So 


U 


1 (red) 


10" 2 


10 48 


1.5 x 10 59 


1.5 x 10 45 


3 x 10 6 


2 (green) 


10" 3 


1Q 48 


1.5 x 10 58 


1.5 x 10 44 


3 x 10 6 


3 (blue) 


10" 3 


10 47 


1.5 x 10 58 


1.5 x 10 43 


3 x 10 7 



magnitude of metal production and so it is only the combi- 
nation of the two, /*, that is of interest. 
The total ionising flux from a halo is 



S = f*M b Q, 



(6) 



where Mb » 0.152 M is the baryonic mass and M is the 
total mass of the halo. We will present results for 3 models, 
listed in tabled Model 1 has the highest ionising flux; model 
2 has a smaller flux but lasts for the same length of time; 
model 3 has an even smaller flux but lasts for longer so that 
the total number of ionising photons produced is the same 
as for model 2. 

Once the stars have formed, the ionising photons will 
begin to evaporate the rest of the halo and make their way 
into the surrounding IGM. Each star cluster produces 



f* M b 
ran 



(7) 



ionising photons, where Mb is the baryonic mass. In the 
absence of recombination, this is sufficient to to ionise the 
hydrogen in a region of baryonic mass 



M bl = N-, 



(8) 



where /xh ~ 1-36 is the mean molecular mass per hydrogen 
nucleus. For the star-formation efficiencies and top-heavy 
initial-mass function that we consider here, there are more 
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than enough photons to ionise any neutral gas within the 
star-cluster: 

^= 80000/„ M h~1.1x10 5 /*- (9) 

We next consider whether it is correct to neglect recom- 
binations. The photon flux required to maintain ionisation 
of the halo (at the mean halo density) is given by 



Sha 



47T 3 2 

— R n n TL, 



(10) 



where R is the radius of the virialised halo, tih is the com- 
bined number density of all spe cies of hydrogen, and 1Z is the 
recombination rate (see, e.g. iHutchings et alJ Eo02*l . The 
value of Shale would be higher if we were to take into ac- 
count clumping of the gas. On the other hand, there are two 
effects that will tend to lower SWo: for high-temperature 
halos not all the gas will be neutral; for low-temperature 
halos the gas will be raised to a temperature that exceeds 
the virial temperature and so will tend to escape from the 
halo — the sound-crossing time for a gas at 10 4 K is of or- 
der 1.0 x 10 7 yr for a 1O 6 M0 halo at an expansion factor 
a = 0.05. We assume that these effects will roughly cancel 
and set the nett photon flux that escapes the halo equal to 
S C sc = S 5haio = fcsc S. Here 



Sh a 



S 



0.12 



(0.O5) 



So 



1.5 x 10 44 phs~ 



where we have set 1Z equal to the recombination rate for a 
10 4 K gas. Of course, / csc is not allowed to drop below zero. 
The number of ionising photons that escape the source halo 
is 



N, 



67, esc 



/esciV, 



fry ■ 



(12) 



Escaping photons are now free to propagate into the 
inter-halo medium (aka inter-galactic medium or IGM) and 
irradiate nearby halos. At the mean density of the IGM, 
the Stromgren radii for any ionised regions are very large, 
and only model 1 produces enough photons to ionise out 
to the Stromgren radius, and then only at very early times, 
a < 0.02. A better picture is that of a bubble of ionised 
gas whose outer radius grows with time until the ionising 
source switches off. To a good approximation then, and for 
simplicity, we assume that recombinations in the IGM are 
negligible. 



5.1.2 Numerical methodology 

The local feedback is implemented as follows. First the tree 
is scanned for all star clusters and a list is generated, in order 
of decreasing star formation redshift. Starting with the first 
cluster, we work up the tree looking at the baryonic mass of 
successive parent halos, M(, iPar until the last halo for which 

M 6 , par < M b + M bl>eac . (13) 

The sub-tree below this parent halo (see Figure |HJ defines 
the extent of the ionised region and the cooling times of 
all halos within it are recalculated taking into account the 
amplitude of the flux and the time for which it acts. 

The tree provides limited information about the spatial 
distribution of halos. However, we know that halos are con- 
fined within a common parent halo and that the parent will 



Merger Tree s . 




First star-cluster 



Figure 8. Schematic view of the merger tree under internal feed- 
back. Once the first star cluster is located, the code calculates 
how many levels up the tree has to go to re-calculate the cooling 
times of each halo under the new root, this time under the influ- 
ence of the radiative flux coming from that first star cluster. In 
this example, the feedback region reaches two levels up the tree. 
Then, the whole sub-tree cooling times will be re-evaluated. 



not have collapsed at this time. So we take the separation of 
the star cluster and each neighbouring halo to be equal to 
the radius of a sphere at the mean density at that time that 
encloses a mass equal to that of the first common parent: 



Rn 



3 M, 



tr.par 



4-7T SlfcO pcQ 



1/3 



(14) 



where a is the scale factor at the time of star formation and 
the other quantities have the usual meanings. Note that the 
value of M„ iP ar in equation ll4l will vary depending upon how 
far one has to travel up the tree to find a common neighbour. 

The flux density, F, at a distance -R par from the source 
is given by: 



F 



AtxRI^ 



(15) 



We need to convert this into an equivalent value of J21 for 
input into our chemical evolution code. To do this we inte- 
grate the spectrum given in equation 2] over all frequencies 
and angles 



F = 4tt 



hv 



du - 



■in 



J21, 



(16) 
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Figure 9. Histograms of the star formation redshifts of primor- 
dial star clusters averaged over a large number of realisations. 
The colours correspond to different models of local feedback as 
listed in table 01 



where hvn = 13.6eV is the energy of H ionisation. Combin- 
ing these two equations allows us to express the ionising flux 
as an equivalent value of J21 for an isotropic radiation field 
(the directionality of the radiation field is unimportant in 
this context). We use the value of J21 and duration of the 
ionisation as inputs to the chemistry code to obtain the new 
cooling time of each neighbouring halo in the sub-tree, using 
the heating and cooling processes explained in Section [2.31 
The list of remaining clusters is then re-ordered using 
the new star-formation redshifts and we continue with the 
next halo in the list. The process is repeated until we reach 
the bottom of the list or until the next halo in the list is not 
able to cool in a Hubble time. 

5.2 Results 

The first stars are widely predicted to form in high density 
regions of space. Consequently, in this section we present 
results of simulations with internal radiative feedback in a 
region for which the root halo corresponds to a positive 3-cr 
density fluctuation. 

Figure [5] plots histograms of star formation redshifts 
for primordial halos for the three cases shown in Table |21 
averaged over a large number of realisations. It is evident 
that the redshift evolution is markedly different for the three 
curves and we shall discuss each in turn. 

First of all, in figure [T01 we compare model 3 with the 
case of zero flux. We can clearly see that the evolution of the 
two curves is identical up to a point, after which model 3 
drops away dramatically. This sudden change can be under- 
stood by examining equation 1111 We are interested in the 
redshift at which the ionising photons first begin to escape 
the halo. If we set / csc equal to zero and solve for the scale 
factor, a, we find that, for this particular model, photons 
do not break out until a redshift z ~ 18, in agreement with 
what is seen in figure [TU1 Also shown are the contributions 
from the two different halo generations, discussed in previ- 
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Figure 10. Histogram of star-formation redshifts for primordial 
haloes. Here we compare model 3 (blue) with the the no flux 
case (grey). Also show are the contributions from Generation 1 
(orange) and Generation 2 (magenta) halos. 



ous sections. The escaping photons have had a devastating 
impact on the surrounding halos, particularly the smaller, 
Generation 1 halos. As a result there is a rise in the number 
of Generation 2 halos because the reduced contamination at 
early times allows more massive halos to cool as primordial 
objects. 

In models 1 & 2 the photons are able to escape the 
halo at much earlier times (before any objects have been 
able to cool). As such, the first primordial objects to form 
will immediately begin to influence their surroundings. This 
explains the much greater reduction in Generation 1 star 
clusters for these models compared to model 3, seen in figure 

m 

Figure tTTl plots the fractional mass per dex of primordial 
objects as a function of temperature for our highest flux 
case (model 1). As expected, the higher flux has suppressed 
the the small, low temperature halos at high redshift, thus 
reducing the amount of early contamination. Once again we 
see that there has been an enhancement in the number of 
high temperature, high mass halos, the distribution of which 
is reminiscent of that seen in figure 2] for the high values of 
J21 in our global model, particularly J21 = 10 -2 . Indeed, 
the average J21 values received by halos in this model are in 
the range J21 = 10 -3 — 10 -2 , consistent with our previous 
results. 

With ten times fewer photons, we expect model 2 to be 
less destructive than model 1 at higher redshifts and figure 
|U] confirms this fact. The green curve shows more primordial 
objects early on which consequently reduces the number of 
Generation 2 objects that form. Interestingly, this model 
shows an approximate balance between the two generations 
with a roughly constant formation rate of primordial halos 
between redshifts z ~ 10 — 22. 

The mass fraction of stars contained in primordial star 
clusters for all the models presented here remains relatively 
constant, varying from 0.06 to 0.13. However, the mass func- 
tions vary substantially as do the star formation histories. 
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Figure 11. Fractional mass per dex of primordial objects as a 
function of temperature, for model 1 (red) and the no flux case 
(grey). 



As the specific ionizing flux increases, the balance moves to- 
wards later star formation and more Generation 2 primor- 
dial star clusters. The positive feedback into Generation 1 
clusters seen in figure 0] for intermediate values of J21 lasts 
for too short a time to be noticable. 



6 DISCUSSION 

This paper looks at the impact of radiative feedback on pri- 
mordial structure formation. This is done in two ways. 

The first part investigates the properties of primordial 
objects under a global UV background. The merger tree is 
illuminated by a constant and isotropic radiation field of four 
different intensities, parametrised by a constant value of J21: 
1(T 10 , 1(T 5 , 10" 2 and 10. It seems more plausible that any 
background radiation field would gradually build up over 
time with the formation of more and more primordial stars. 
Consequently, we also investigate a time-dependent build- 
up of_th*_backgjwnrl^_ux using an extension of the model 
of iKitavama et alJ i200Ct) . The section of the paper uses a 
mean-density region of space as the background radiation 
field is assumed to come from external sources within higher- 
density regions. 

The effect of a constant UV field on the halo popula- 
tion is not a trivial one as both positive or negative feedback 
can arise from different choices of the flux amplitude J21. 
The cooling of a low-temperature primordial gas is almost 
completely dominated by the release of energy from roto- 
vibrational line excitation of H2. But if a radiation field is 
present, H2 is easily destroyed by Lyman- Werner photons 
(11.2eV< hv < 13.6eV). On the other hand, the formation 
of H2 can be enhanced by an increase in the ionization frac- 
tion produced by a weak ionising flux, as electrons act as a 
catalyst for the formation path of H2. At modest flux levels 
of J21 ~ 10 -5 the nett effect is to boost cooling in the first 
star clusters. A similar result has previously been found by 



lHaiman et alJ Jl996ft. iRicotti et alJ <200ll) . IKitavama et all 
fcOOll) and lYoshida et alJ i2003t) . 

Negative feedback is produced not only by photo- 
dissociation of H2; at high flux levels the dominant effect 
comes from heating associated with photo-ionisation of H. 
For fluxes of J21 > 10 -2 we find that molecular cooling 
is ineffective and only haloes with virial temperatures of 
T v i r > 14000K are able to cool in a Hubble time. Because of 
reduced contamination from star-formation in low-mass ha- 
los, strong radiation fields can increase the number of high- 
mass primordial star clusters. 

The second part of this paper dealt with a model in 
which the radiative feedback is localised. That is, star clus- 
ters irradiate their surrounding area, changing the cooling 
properties of those primordial objects that are inside their 
ionisation spheres. For this model we considered only a high 
density region corresponding to a 3-ct fluctuation since the 
first objects are thought to form in regions of high overden- 
sity. 

The effect of the feedback depends mainly upon the 
specific ionising flux, averaged over the mass of the halo. 
When this is low, most or all of the photons will be used 
up in maintaining the ionisation of the halo. For the par- 
ticular model that we describe in this paper, equation 1111 
relates the escape fraction, / osc , of ionising photons to the 
specific ionising flux, So, and redshift. A value of So below 
10 43 phs~ 1 Mq^ 1 will reduce / osc to zero until the redshift 
drops to about 16, corresponding to the peak in the produc- 
tion rate of primordial halos per unit redshift. Although the 
precise number will be model dependent, we regard this as 
a fiducial value below which feedback will be ineffective. 

Higher values of So result in shift from primordial star 
clusters away from Generation 1 (low virial temperature) 
towards Generation 2 (high virial temperature). Unlike the 
case of a global ionisation field, Generation 1 clusters are 
not erradicated completely, because some must form in order 
to provide the feedback. However, a specific ionisation flux 
of So = 10 45 phs^Mo" 1 is enough to swing the balance 
strongly in favour of Generation 2 star clusters. 

This paper makes a number of advances on our previous 
modelling, most notably the use of an improved merger tree 
that does not restrict halos to factors of two in mass, and 
the introduction of a radiation field. The former results in a 
reduced mass fraction of stars in primordial halos with the 
bias shifting more in favour of Generation 1; however, the 
latter moves the bias back the other way. The conclusions 
of HSTC02 and ST03 remain valid in that there could be a 
substantial population of primordial star clusters that form 
in high-mass halos dominated via electronic cooling. 

Further improvements to our model are possible. Al- 
though we do not expect these to change our qualitative 
conclusions, they will be important for making precise quan- 
titative predictions about the number density and compo- 
sition of the first star clusters. We mention some of them 
below. 

We assume in this paper that the internal structure of 
halos is unaltered between major merger events. However, 
it is possible for halos to increase their mass substantially 
through a succession of minor accretion events, and this will 
release gravitational potential energy and lead to heating 
jYoshida et alJl2003t iReed et al.ll2005l) . We intend to incor- 
porate this in future work. 
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No attempt has been made to distinguish between 
star-formation and feedback in primordial star clusters and 
later generations. It is expected that the initial mass func- 
tion should alter c onsiderably once the metallicity reaches 
abou t 10~ 3 5 Zrr, jBromm fc Loebl 120031: ISchneider et alJ 
120021: ISantoro fc ShuUll2005l) and that the spectrum of the 
ionization field would become softer and its normalisation 
lower. 

We have neglected the effect of stellar winds and su- 
pernovae, both of which help to heat the surrounding gas 
and pollute it with metals. The effect of winds is fairly lo- 
calised, but supernovae can lead to superwinds that expel 
enriched material from the star cluster. It is possible that 
this could affect neighbouring halos; more likely it will sim- 
ply ensure thorough mixing of metals throughout a common 
parent halo. Numerical simulations will be required in order 
to model this process with any degree of realism. 

Finally, it would be valuable as computing resources im- 
prove to simulate larger regions of space that would provide 
a more respresentative section of the Universe. This would 
allow us to distinguish between the early history of galaxies 
that are located within clusters and in mean density regions 
of space, for example. It would also allow for both localised 
feedback and a self-consistent build-up of the background, 
global radiation field. 
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The appendix below describes how to go from cross- 
sections to reaction rates and heating terms. 



APPENDIX A: PHOTO-IONISATION AND 
PHOTO-DISSOCIATION INTEGRALS 

The rate at which photo-ionisation or photo-dissociation re- 
actions occur is given by: 

h=4-K / a u ,i ^ du, (Al) 

where /if t h is the threshold energy for which photo- 
ionisation (or photo-dissociation) is possible and a v ^ is the 
frequency dependent cross-section of the i th reaction. The 
cross-section and the threshold energies were taken from Ta- 
bic J„ is the specific intensity given in Equation 2] 



APPENDIX B: HEATING TERMS 

The energy per particle that a photon of energy hl> trans- 
fered to an electron in an atom or ion is hv — huth- Therefore, 
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the energy per second per unit volume transferred to the gas 
(heating terms or heating functions) is: 

A hc at,i = 4-7rni / o v ,i J v — dp, (Bl) 

•'"th 

where m is the number density of the dissociated species. 

The heating of a primordial gas comes primarily from 
photo- ionisation of H, He and He + , but there is also a small 
contribution from photo-dissociation of H2. 
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